Spin wave dispersion softening in the 
ferromagnetic Kondo lattice model for manganites 
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Spin dynamics is calculated in the ferromagnetic (FM) state of the generalized Kondo lattice model 
taking into account strong on-site correlations between e g electrons and antiferromagnetic (AFM) 
exchange among ti g spins. Our study suggests that competing FM double-exchange and AFM 
super-exchange interaction lead to a rather nontrivial spin-wave spectrum. While spin excitations 
have a conventional Dq 2 spectrum in the long-wavelength limit, there is a strong deviation from the 
spin-wave spectrum of the isotropic Heisenberg model close to the zone boundary. The relevance of 
our results to the experimental data are discussed. 



The revival in the study of manganites has led to ex- 
perimental re-examination of their different properties 
and interplay between nontrivial magnetism and highly 
anisotropic charge transport. One of the puzzling fea- 
tures is the non-universality of the magnetic and trans.- 
port properties of FM manganites at doping x ~ 0.3a. 
According to the conventional theory of the Double Ex- 
change (DE), the spin dynamics of the FM state that 
evolves at temperatures below the Curie temperature Tc 
is expected to be of nearest-neighbor Heisenberg type. 
This picture seems to be indeed reasonably accurate 
for manganites with high value of Tea Recently, how- 
ever, several experimental results have shown a strong 
deviation of the spin-wave dispersion (SWD) from the 
typical Heisenberg behavior. The unexpected soften- 
ing of the SWD at the zone boundary has been ob- 
served in several manganites with Iova but very different 
T c : Pro.esSfio.srMnOa (T c =301 K)B, Lao.7Cae. 3 Mn0 3 
(T c =242 K)d and NdcrSr^MnOa (T c =198 K)B. These 
observations are very important as they indicate that 
some aspects of spin dynamics in manganites have not 
been entirely understood yet. During last few years this 
problem was discussed by several theoreticians. Khali- 
ullin and Kilian in Ref . ^| proposed a theory of anomalous 
softening in FM manganites based on the modulation 
of magnetic exchange bonds by orbital degree of free- 
dom of double-degenerate e g electrons. They found out 
that charge and coupled orbital-lattice fluctuations can 
be considered as the main origin of the softening phenom- 
ena. Soloviev and Terakura in Ref. fj] have argued that 
the softening of spin wave at the zone boundary and the 
increase of the spin-wave stiffness constant with doping 
have purely magnetic origin. 

The DE mechanism as an_explanation of ferromagnetic 
order in doped manganitestj was propos/ad shortly after 
the first experiments in these compounds! In this model 
it is assumed that the hopping of e g electrons between 
neighboring sites is easier if the local spin on the sites 
are parallel, so an effective ferromagnetic coupling be- 
tween the local spins is induced by the conduction elec- 
trons lowering their kinetic energy. For quite a long time 



it was believed that DE model gives full understanding 
of the FM metallic state in manganites. Later it was 
shown both from theory and experiment that, in fact, it 
provides only a qualitative explanations and it is neces- 
sary to extend DE model by taking into account the or- 
bital degree of freedom of double-degenerate e g electrons, 
phonons, AFM superexchange interaction between local- 
ized spins , etc. Particularly, for the Curie temperature 
the value estimated by pure DE was much larger than 
the experimentally observed one. Taking into account 
AFM superexchange interaction among localized spins 
of t2 g electrons gives possibility tOj-diminish the value of 
Tc in the framework of DE modelEl Recently, a Monte 
Carlo study of the Ferromagnetic Kondo Lattice Model 
(FKLM) for doped manganites (see, Ref. |ll|) confirmed 
the scaling of the Tc suppression with a bandwidwidth 
narrowing induced tw the antiferromagnetic frustration. 
In our recent papertHl we have also shown that AFM su- 
perexchange interaction plays very important role in de- 
forming the phase diagram in manganites. 

The aim of our work is to show that the competition 
between AFM superexchange interaction among local- 
ized spins and double-exchange interaction among itin- 
erant electrons is one of the mechanisms responsible for 
the spin wave softening phenomena. 

In the present paper we continue to study the SPfn dy- 
namics within the generalized one-or bitaLFMKLB. The 
effective Hamiltonian of FMKL modelEj'ta can be written 



H = ~ H$ij)c ] (i)c(j) + ™t(*W0 ; ) 



J ff ^S(i)s(i) + l/2^J«S(i)S(i) , 



(1) 



where c(i) and c'(i) are annihilation and creation opera- 
tors for electrons in the spinor notation at site R(i); 



c T (i) 
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\x is the chemical potential; n a (i) = c^(i)c a (i) is the den- 
sity operator of electrons with spin cr; the spin operator 
is given by s&(i) = ^c^ (i)<Jkc(i), (k = 1,2,3) are the 
Pauli matrices; S(i) is the localized Mn core spin of t^g 
electrons with S — 3/2; tij is the nearest-neighbor hop- 
ping parameter; Jh is the FM Hund coupling interaction 
among the localized and delocalized spin subsystems, Jij 
is the nearest-neighbor antiferromagnetic superexchange 
interaction among localized spins. 

The magnetic properties of the above given Hamil- 
tonian has a dual character: there exists a purely 
Heisenberg-like contribution from the exchange interac- 
tion among localized spins and the itinerant contribution 
driven by electronic excitations. The itinerant contribu- 
tion leads to the effective FM exchange interaction be- 
tween localized spins due to the DE mechanism. To treat 
appropriately this dual character of the magnetism one 
should describe both contributions on the equal level of 
approximation as well as retain SU(2) spin symmetry of 
the Hamiltonian. To this end, we enmloy the Compos- 
ite operator method (COM) approacblij by introducing a 
new set of operators which well describe collective excita- 
tions in the system and are constructed from the original 
electron operators. These composite operators can have 
fermionic or bosonic character. They are created by in- 
teractions among the electrons and the localized spins, 
and, therefore, their properties will be determined by 
the dynamics and boundary conditions and must be com- 



puted self-consistently. 

First we start by discussing the electronic excitations of 
the system, responsible for the itinerant contribution to 
spin dynamics. We introduce the 4-component fermionic 
basis: 



(2) 



where and n a {i) are Hubbard composite excitations: 

(3) 



€a(i) = c CT (i)(l - n s (i)) 
f] a (i) = c a (i)n 9 (i) . 



The r] excitation describes an electron restricted to move 
on sites already occupied with an electron of opposite 
spin whereas £ demands that there be no prior occupancy 
on the site. 

Let us consider the equation of motion for the field 



(4) 



where is the current operator. By considering the 
Hamiltonian given by (Q) , we obtain the following ex- 
pression for the current operator: 



( - 6tcf(i) - 6^ T « - J H /2 Ki(i)S-(<) + ^(i)S z (i)] \ 

-QjL-U) m (i)+Gt7ri(i)-J H /2 [Vi(i)S-{t)+ m (i)S'(i)] 
-/*;(*) - 6tc?(i) - 6t n (i) - J H /2 [$(i)S+ (i) - 
V -([1- 17)^(1) +6t n (i)- J H /2 [m(i)S+{i)-Vim z (i)} J 



(5) 



where n(i) = i<T M n At (i)c 0! (i) + c(i)c aj( c(i) is a higher or- 
der composite field; n^(i) = c^(i)u^c(i) - the number 
(fx = 0) and spin (/i = 1,2,3) density operator with the 
notation = (l,er), er M = (— 1,<t), a being the Pauli 
matrices; and the notation c a (i) = ■ a(i,j)c(J) stands 
to indicate the field c on the first neighbor sites with 
F.T.a(i,j) = ^(cosk x + cosk y + cosfc 2 ) (abbreviation 
F.T. stands for Fourier transform). 

Formally, we can write the nth element of the current 

as 

3n(i) = ^2 ^2 £nm(hj)i>m(j) + S 3n(i) ■ (6) 
j m 

By projecting out the correction part 5j n (i), we linearize 
the Heisenberg equation of motion as 

d 

l g:^n(i) = 22£nm(hj)lpm(j) ) ( 7 ) 
3 

where the energy matrix s(i,j) is calculated by means of 
the equation 



e{i,3)=Y, m ( i ' l ) I ~ l ^3) ■ ( 8 ) 
i 

To shorthand the notation we defined the normalization 
matrix and the frequency matrix m(i,j): 

=< {il>{i),^{j)} E .T. > , 

=< {j{i),ifi{j)}E.T. > , 

where E.T. indicates the anticommutator evaluated at 
equal-time; the bracket < ... > denotes the thermal av- 
erage on the grand canonical ensemble. 

Straightforward calculations give us expressions for the 
Fourier transforms of the matrices / and m. The normal- 
ization matrix is diagonal and has the following matrix 
elements: 

In = 1- < nx >= 1 - n/2+ < s z > , 
J22 =< nx >= n/2- < s z > , 
hs = 1- < n T >= 1 - n/2- < s z > , 
I44 =< n\ >= n/2+ < s z > ; 
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where n =< n > is the particle density and < s z >= 
1/2 < n 3 >. The m matrix is hermitian and the ele- 
ments different from zero are given by 

mn = -nhi - 6iA x - 6ta(k)(Jn - I 22 +pt) 

- Jff(a + 6 2 ) , 

mi2 = 6iA^ +6ta(k)(-/ 2 2 +Pt) > 

"^22 = -(m _ ^0-^22 - 6tAj - 6ia(k)p T + J H b 2 , 

^33 = -M-^33 - 6tA T - 6ia(k)(/ 33 - 7 44 

- Jff(-a + 61) , 

m 34 = 6iA T + 6ia(k)(-J 44 , 

m 44 = — (fj, — U)In — 6iA| — Qta(k)pi + Jnbi , 

(10) 



composite fields. If these composite fields belong to the 
fcrmionic or bosonic basis they can be expressed in terms 
of the corresponding GF. However, it may happen that 
some of the parameters are expressed as expectation val- 
ues of higher-order composite fields that do not belong to 
the basic set and we need to evaluate by imposing some 
symmetry requirements, as it will be discussed later. 

As we have already mentioned the spin dynamics of the 
present model is governed by both localized and itinerant 
spin subsystems. To treat these two contributions at the 
same level of the approximation we consider a bosonic 
sector by defining a basis of composite fields as 



B(i) 



s+(i) 
S+(i) 



(14) 



where a =\ (1) or c =J. (2). Matrices m(k) and I involve 
correlation functions which we defined as 

A<r = < > - < > , 

Pa = j[< n£(*K(») > +2(-l) ff < n a (i)m(i) >] 

-<c 1 {i)c 2 ^[c\{i)c\{i)r> , (11) 
a = 1/2 < S z > , 

K = l/2(< s-(i)S+(i) > -(-I)* < n a (i)S z (i) >) . 

The parameters A CT and p a are static inter-site corre- 
lation functions which describe, respectively a constant 
shift of the bands and a bandwidth rcnormalization in 
itinerant subsystem. 

The retarded Green function (GF) can be written as 



The procedure then is similar to the fermionic sector; the 
propagator for the boson field B(i) will contain another 
set of parameters which must be self-consistently deter- 
mined. After some algebra the following expressions of 
the bosonic normalization matrix 

J(k) =F.T. < [B(i),B\j)] E . T > 

and the frequency matrix which is given by the commu- 
tator of the current and the bosonic fields 

ro(k) - F.T. < [i^B[i),B^{j)} > E .T. 
can be easily derived: 



s F (i,j) =< W(*),v t 0')}>= 

. ; j d 3 k J uexpiHn-r^-Mti-ti) S F (k,Lu) , (12) 



in 



where the k-integration is over the Brillouin zone (BZ) 
and Q is the inverse volume of the BZ. By means of the 
linearized Heisenberg equation (7) the Green's function 
has the following expression: 



S F (k,Lj) = 



1 



w-m(k)7" 1 (k) 



J(k) 



or, in the spectral form: 



S F (k,L0) 



4 

E 



<7*(k) 



uj-Ei (k) + iS 



(13) 



where the energy spectra Ei (k) and the spectral functions 
CT*(k) can be easily evaluated. The energy spectra Ei(k) 
are the eigenvalues of the matrix e(k). Since they depend 
on a set of external parameters such as electron density 
n, temperature T, Coulomb interaction U, Hund cou- 
pling Jh and AFM Heisenberg exchange interaction J, 
and a set of internal parameters (/1, m, A a ,p a , a, b\, b 2 ), 
they must be calculated in a self-consistent way. Inter- 
nal parameters are expressed as expectation values of 



In = 2 < s z > , 

/12 = I21 = , (15) 
I22 = 2 < S z > ; 

^11 = Jnd + 6tto(l — a(k)) , 

mi2 = m 2 i = - Jnd , (16) 

m 2 2 = Jnd - 6J/(1 - a(k)) , 

with coefficients defined as 
w =< c*{i)c a {i) > , 

d = 2 < s z {i)S z (i) > + < s+(i)S-(i) > , (17) 
/ = 2 < S°{%)S x {i) > + < (S+) a (i)S-(i) > . 

It is important to emphasize that w is the kinetic energy 
and, hence, the effective FM exchange scales with the 
kinetic energy. 

The retarded Green's function for the boson field B(i) 
can be expressed as S B (i,j) =< R[B{i)B^{j)\ >. In the 
polar approximation the bosonic GF can be written as 



cr*(k) 



f w-25j(k) + i<J ' 



(18) 



where o % (k) and Ei (k) are bosonic spectral functions and 
the energy spectra, respectively. Ei(k) are two branches 
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of spin excitations: acoustic and optical modes. In the 
long-wavelenght limit the accoustic mode is -E'i(k) ~ 
Dk 2 , and the spin stiffness scales as D ~ (ztw — J) which 
clearly shows the competitions of the effective, induced 
by the electron hopping, FM exchange Jfm and AFM 
superexchange J. 

There are new internal parameters appearing in the 
bosonic sector. From the analysis given above, the 
bosonic and fermionic GF depend on 12 internal param- 
eters: 



< s z >, <S Z >, < s+(i)S~{i) >, fj, 
A CT , p a , w, f, < n a (i)S z (i) > . 



(19) 



Now we need to construct a closed self-consistent 
scheme for fixing these parameters. The static correla- 
tion functions (CF) Ci rn (i,j) can be computed from the 
corresponding retarded GF by the fluctuation-disipation 
theorem: 



C F(B) (M) 
" -to 



( 



< B(i)B^{j) > > 
dwe ik(ri ~ rj) "^ (ts 



(20) 



-to) 



C(k,w) 



where 



(I + tanh(cj/2T))ImS* f (k, uf) 
(I + coth(w/2T))ImS s (k, lu). 



(21) 



For convenience, we will use the following notation for 
the static CF: C lm = C lm {ify Cf m = C lm (i,i a ). 

It should be pointed ouO that one can use Eq.(|2l|) 
for calculation of correlation functions C(k,u>) in terms 
of the commutator bosonic GF only for ergodic systems 
when there is no 5(u>) term in the spectral density. Other- 
wise one has to use the anticommutator or causal GF for 
bosonic operators to take into account nonergodic terms. 
In our case of FM ordered system with < S z ^ > 
and < s z =/= > the bosonic GF ( |l8| ) describes the one- 
particle, spin-wave-like excitations which are ergodic for 
all q ^ 0. 

The first self-consistent equation is given by the defi- 
nition of the magnetization per site: 



The kinetic energy w can also be expressed through the 
linear combination of GF matrix elements: 



-( E + E )s 

n , m — 1 n.m— 3 



Fa 
nm 



(25) 



The other parameters are not strictly bound to the 
dynamics and they must be evaluated by other means. 
This aspect, is a general property of the Green's function 
formalismta. The equations of motion are not sufficient 
to completely determine the Green's functions which re- 
fer to a specific choice of the Hilbert space. The infor- 
mation concerning the representation must be supplied. 
As discussed in Refs. [l7| and |l^, this information is sup- 
plied by determining the parameters not bounded by the 
dynamics in such a way that the Hilbert space has the 
right properties to conserve the relations among matrix 
elements imposed by symmetry conditions. In particu- 
lar, the main attention should be paid to the fact that 
the chosen represention does not violate the symmetry 
required by Pauli principle. Let us note that by Pauli 
principle we mean all relations among considered oper- 
ators dictated by their algebra. The Pauli principle for 
the electron propagator requires that 



°12 
fiF 

n F _ 
O n - 



= 
= 
°33 



(26) 



There is one more self-consistent equation which connects 
matrix elements of fermionic and boson GF: 



fit} 



B _ r <F 



(27) 



To close our self-consistent scheme, we have to do 
some approximation and we use the following decoupling 
scheme: 



< n a (i)S z (i) >~< n a {i) >< S z (i) > 



(28) 



f 



(n T 



f 



la 



(22) 



The boundary condition which fixes the particle number 
as an external parameter gives us the equation for fixing 
the chemical potential /i: 



2 - C£ 



nF 

L/99 



/-i t _ n 
°33 °< 



F 

44 • 



(23) 



The other three parameters are expressed through the 
matrix elements of the CF ([H]) as follows: 



a t = c£ a - a 



22 

A r*<Fot nFa 

— u 33 — °44 

<*+(*)£-(*)>= C? 2 



(24) 



Having solved the system of the self-consistent equa- 
tions we can study various properties of the FKLM. In 
this paper we will look at the spin wave spectrum. The 
study of other properties will be published elsewhere. 

In Fig. I we plot the dispersion of magnetic excitations 
along various direction in the Brillouin zone, and we com- 
pare it with nearest-neighbor Heisenberg form with the 
spin stiffness D = lim g ^ E(q)/q 2 , where E(q) is the 
acoustic mode of the bosonic spectrum calculated in the 
mean field approximation (see, formula (24) in Ref. [l3]). 
We note that in this approach the acoustic spin mode 
already occurs in the mean-field approximation as the 
result of the full rotational symmetry of the spin system. 
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FIG. 1. Spin wave dispersion obtained in FKLM along 
(tt, 0,0), (7r,7r,0), (7r,7r,7r) directions (solid line). The fol- 
lowing external parameters are chosen: Jh = Wt, J — O.lt, 
T — 0.3i, U = At. Dotted line represent ideal Heisenberg 
spectrum with stiffness D. 



of exchange bonds. If orbitals fluctuate slower then spins 
then the anisotropy imposed upon the magnetic exchange 
bonds by the orbital degree of freedom forms spin dy- 
namics. The presence of Jahn- Teller distortions quenches 
the dynamics of orbitals and can lead to the long range 
orbitally ordered state which corresponds to the mini- 
mal frequency of orbital fluctuations. Then both orbital 
fluctuations and the competition between the double- 
exchange interaction and the AFM superexchange inter- 
action are responsible for the strong renormalization of 
the magnetic excitation spectrum. 

In Pro.63Sro.37Mn03 there is a tendency to the planar 
(x 2 — y 2 ) ferrotype orbital ordering. That means that 
orbital fluctuations mostly renormalize exchange bonds 
in two dimensions. Let us note that magnons in (tt, tt, 7t) 
direction are sensible to all three spatial directions of 
the exchange bonds and thus, their dispersion remains 
unaffected by orbital fluctuations and softening of the 
spectrum in this direction is purely due to the AFM spin 
fluctuations. In our approach we attribute softening phe- 
nomena to the effect of the suppression of the FM order- 
ing only by the AFM superexchange interaction and thus 
on the direction (tt, tt, tt) we estimate ordering correctly. 
In other direction we have to add influence of orbital de- 
gree of freedom to fit experimental results. 

In summary, we have shown, that the spin wave spec- 
trum obtained by the COM analysis of FKLM shows the 
softening close to the magnetic zone boundary. In our ap- 
proach, we have accounted for both the intersite AFM ex- 
change among localized i2g-electron spins and the strong 
intra-atomic Hund coupling among the ti g and e g elec- 
trons and proved that the competition between these two 
interactions can lead to the strong renormalization of the 
magnon spectrum and suppression of the FM ordering. 

We thank G. Khaliullin and G. Jackeli for valuable dis- 
cussions. N.B.P. wishes to acknowledge partial support 
by Unita I.N.F.M. di Salerno. Financial support by IN- 
TAS Program, Grants No 97-90963 and No 97-11066, are 
acknowledged by N.B.P. and N.M.P.. 



At small momenta the spectrum exhibits conventional 
Heisenberg behavior. With increasing momentum the 
spectrum shows the softening which has the maximum 
value close to the (tt, tt, tt) AFM vector. However, in ex- 
periment the maximum softening was observed in (tt, 0, 0) 
and (tt,tt,0) directions. In order to reproduce the ex- 
perimental data it is thus important to consider orbital 
degree of freedom as it was considered by Khaliullin and 
Kilian (see Ref. ||). How strong the orbital fluctuations 
modulate exchange bonds between Mn ions depends on 
their characteristic time scale. If the typical frequency 
of orbital fluctuations is higher than the one of spin fluc- 
tuations, the magnon spectrum is renormalized only due 
to AFM fluctuations and then the orbital state effects 
the spin dynamics only by restoring the cubic symmetry 
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